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ABSTRACT 


Exact  solutions  for  the  transient  temperature  distribution  and 
the  stored  energy  in  an  infinite  plate  of  finite  thiclmess  are  presented 
for  the  case  of  different  convective  environments  at  each  face  of  the 
plate.  The  solution  is  general  and  contains  numerous  limiting  cases, 
Including  that  of  steady  state.  Eigenvalues  axe  given  for  many  . 
combinations  of  the  system  Biot  numbers  for  the  initial  response  period. 
An  example  is  presented  to  illustrate  the  application  of  the  solution 
to  the  practical  problem  of  a  rocket  engine  diffuser. 
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SECTION  I 
INTRODUCTION 


The  transient  response  of  an  infinite  flat  plate  of  finite 
thickness  has  been  analyzed,  for  many  cases  (Refs.  1  through  5).  However, 
the  most  general  solution  for  convective  environments  (different  surface 
heat-transfer  coefficients  and  thermal  environments  on  each  side  of  the 
plate)  is  not  available,  although  the  possibility  of  the  solution  is 
mentioned  in  Ref.  1.  The  transient  response  of  a  plate  subjected  to 
unsyrametrlcal  boundary  conditions  is  very  important  in  many  analyses. 

For  example,  the  transient  time  is  the  prime  period  of  interest  in 
evaluating  the  behavior  and  application  of  structures  subjected  to 
unsymmetrical  boundary  conditions,  such  as  exhaust  gas  diffusers  for 
simulating  the  hi^  altitude  environment  of  rocket  engines,  rocket  engine 
nozzles,  ejectors,  tunnel  walls  of  high  temperature  short  run  time  test 
facilities,  nozzles  of  intermittently  operated  rockets,  and  components  of 
aircraft  and  missiles  in  high  speed  flight.  In  many  of  the  above  cases, 
the  engineer  Is  ultimately  interested  in  predicting  coolant  flow  rates 
necessary  to  keep  the  wall  within  structural  and  material  temperature 
limits.  It  is  reasonable  to  expect  that  the  coolant  rate  necessary 
for  a  short  time  test  or  exposure  may  be  of  a  reasonable  magnitude, 
whereas  the  coolant  rate  necessary  for  steady-state  operation  may  be 
completely  unreasonable  in  some  of  the  above  applications. 

This  analysis  presents  an  exact  solution  for  the  temperature 
response  in  a  solid  bounded  by  two  parallel  planes  with  unsyrametrlcal 
boundary  conditions .  Implicit  in  the  solution  is  the  capability  of 
predicting  a  coolant  flow  rate  necessary  to  keep  an  exposed  wall  within 
structural  and  temperature  limits. 

SECTION  II 

ANALYSIS 


2.1  PHYSICAL  SYSTEM 

The  physical  system  considered  in  this  analysis  is  shown  in  Fig.  1, 
An  infinite  plate  of  finite  thickness  Z  Is.  initially  at  a  uniform 
temperature  f(x)  throughout.  At  time  t  ^  o,  the  face  at  x  =  0  is  exposed 
to  a  high  temperature  convective  environment  at  temperature  v  .  Similarly, 

for  time  t  >  o,  the  face  at  x  =  i  is  exposed  to  a  lower  temperature 

convective  environment  at  temperatiire  v  .  Assume  that  the  surface  heat- 

2 

transfer  coefficients  h  and  h  are  uniform  and  constant  at  x  =  o  and 

1  2 

X  =  i,  respectively.  The  thermal  conductivity  and  thermal  diffusivity 
are  given  by  k  and  K,  respectively,  and  are  assumed  to  be  independent 
of  temperature  and  position. 


1 


AEDC-TR-66.166 


2.2  MATHEMATICAL  MODEL 

A  'basic  energy  'balance  on  the  plate  shows  that  the  partial . 
differential  equation  describing  the  temperatxure  distribution  in  the 

plate  is  given  by  ^  _  k  ^  ( 1 

This  equation  is  subject  to  the  following  boundary  conditions; 
k|^-h(v-v)=0  atx  =  0 


k—  +h(v-v)=0  atx=i 
Sx  ^  2 

and  the  initial  condition: 


V  =  f(x)  at  t  —  0 


The  above  system  of  equations  can  be  solved  by  many  different 
techniques.  However,  the  principle  of  superposition  is  especially 
convenient  for  this  problem.  Assuming  that  the  solution  can  be  expressed 
as 

v(x,t)  =  u(x)  +  w(x,t)  (2) 

where  u(x)  is  the  steady-state  contribution  to  temperature  and  w(x,t) 
is  the  transient  contribution,  then  u(x)  must  satisfy  the  differential 
equation 

ifs.o  (3) 

dx^ 

subject  to  the  following  boundary  conditions: 

k|ii-h(u-v)=0  atx  =  0 
dx  i'  i' 

k^  +  li(u-v)  =  0  atx  =  i 
dx  2  s 

The  function  w(x,t)  must  then  satisfy  the  partial  differential  equation 

^  =  O^xSi  (l^) 

9t 
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subject  to  the  following  boundary  and  initial  donditions 


Sw 

k  —  -  hv=0  atx*0 
Sx  ^  * 


dw 

k  —  +  hw  =  0  atx=i 
3x  2 


w  =  f(x)  -  u  at  t  ^  0 


It  may  be  readily  shown  that  the  solution  to  the  system  of  Eq.  (3)  1b 


HH  (v  -v)x+Hv  (l+H^)+Hv 
1  11  2  2  2 

K  +  H  (1  +  H  i) 

12  1 


(5) 


h  h 

where  H  s  ^  and  H  = 

1  k  2  k 


It  may  be  shown  that  by  using  the  product-type  solution,  the 
solution  to  the  system  of  Eq.  (h)  is 


n  =  1 


CO 


where 


Z  (x) 

v»  '  ' 


[2(p^  +  H^)j^|p^cos  (p^x)  +  Bln  (P^x)] 

{K  +  ^  *  %]  +  \K  *  KV 


(T) 


where  are  the  positive  roots  of 


(p®  -  H  H  )  6in  (p  =  p  (H  +  H  )  COB  (p^^) 
di2  «  nig  n 


(8) 


Therefore,  the  solution  to  Eq.  (l),  using  the  assumption  of  Eq»  (2), 
is  the  s\jm  of  Eqs.  (5)  S'Od  (6),  or' 


HH(v  -v)x+Hv  (l+Hi)+Hv 
v(x,t)  =  i ^ 

H  +  H  (1  +  H  £) 

IS  1 


(2) 


-K  p®t  I 
'  *<1 


^  Z^(x)e  Z^(x')  [f(xO  -  u(x')]  ax' 


n  ^  1 
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It  is  beyond  the  scope  of  this  work  to  prove  that  Eq.  (2)  represents 
the  unique  solution  to  the  system  of  Eq.  (l)  and  that  Eq.  (2)  is  a  uniformly 
convergent  solution;  uniqueness  and  viniform  convergence  may  he  shown  readily. 


For  simplicity,  let  the  general  solution  (Eq.  [2])  he  modified  hy 
assuming  that 

f ( x)  =  f  ( x’)  =  V  =  initial  plate  temperature  (  9 ) 


Substitute  Eq.  (9)  into  Eq.  (2)  and  integrate  to  obtain 


(v  -  V  )  Th  +  H  H  X 
2  1  L  2  12. 


V  -  V  =  - ^ 

^  H  +  H  (1  +  H  i) 
12  1 


p 

“  (p2  +  h2)  [p  cos  (P  x)  +  H  sin  (P  x)]  H  (V  -  V  )  +  (H  +  H  H  i)(V  -v) 

-f.  2  n  2L^  n  1  nj  1  1  3  12  2 

n  =  1  (p^  +  H^)  ri(p2  +  H®)  +  H  1  +  H  (p^  +  H®)  H  +  H  (l  +  H  i) 

n  L  s'  s\  1  s  L\  i  a  i 


H  (v  -  V  ) 

12  2  1 _ 

K  [\  +  *  H/)] 


H  (V  -  V  ) 
+  -i  ^ 
Pn 


H  H  H  i(v  -  V)  +  (H  +  H  )(v  -  V) 

sin(p  i)  +  ^  - - - - - - - - - 

^n  H+H(l  +  Hi) 


cos(p^i) 


Equation  (l 0 )  will  be  more  convenient  to  work  with  in  a  dimensionless 
form.  Using  the  dimensionless  temperature  ratios  9,  6,  ilf,  based  on 

h  ,  Nt,  based  on  h  ,  and  the  dimensionless  parameter  oi  and  e  ,  Eq.  (lO) 

1  -02  2 

may  be  written 
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The  equation  for  the  eigenvaluea^  Eq.  (8)  In  dimensionless  fonn  becomes 

(12) 


tan  (£^)  =  - 


Equations  (ll)  and  (l2)  are  sufficient  to  determine  the  dimensionless 
temperatiire  distribution  in  an  infinite  plate  of  finite  thicloiess  expjosed 
to  unsymmetrical  boundary  conditions. 

In  addition  to  checking  Eq.  (ll)  for  uniqueness  and  uniform 
convergence,  one  may  also  shov  that  it  reduces  properly  to  various 
"special  cases".  The  Helsler  or  Groeber -type  solution  (face  at  x  =  0 
insulated  or  IL  =  0,  and  N„  finite)  available  in  most  textbooks  on  heat 

transfer  is  readily  obtained  by  letting  =  0  in  Eq.  (ll).  Other 

Pi 

cases,  such  as  both  faces  insulated  (E^  =  =  O),  zero  thermal 

resistance  at  x  =  0  (N^  =  “)j  and  finite  (or  vice  versa),  and 

zero  thermal  resistance  at  both  faces  (Kg  “  ere  all  readily 

obtained  by  proper  reduction  of  Eq.  (ll).  For  the  steady- state  case, 

Eq.  (ll)  reduces  to 


e  (x) 


1  +  Ng^(x/i) 

1  +  a  + 


(13) 


The  total  energy  stored  in  the  plate  per  unit  area  Q,  in  time  t, 
is  given  by  ^ 


Q  =  k(T 


J  (dt  X  =  i 

3x 

•x  =  0 

dt 


The  maximum  energy  stored  in  the  plate  per  xinit  area  is  defined  as 

=  pc  i  (v^  -  v^)  (15) 

Substituting  Eq.  (ll)  into  Eq.  (1^4-),  dividing  by  Eq.  (I5),  and 
performing  the  Indicated  operations  gives  the  ratio  of  the  total  heat 
flow  into  or  out  of  the  plate  in  time  t,  to  the  maximum  energy  of  the 
plate,  as 


±  _  V  l—\ 

%  1  +  «  +  jLj  \e^) 

n  =  1 

|i^  1“  'I'  +  (1  +  sj  -  sin  (e^^)  -  d  +  (I  +  a)i|  cos  (l+a  +  Ng^) 

(16) 


cos  (e^) 

-  Bln  ( e 

I'-l 

i.&t 

V  6  / 
\  n  / 

1-^ 

•&S] 

[fe 

r 

a  ijr  + 

^2  Kt 

1  -  e' 
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Equations  (ll),  (12),  and  (l6)  are  sufficient  to  determine  the 
temperature -time  history  and  the  total  heat  flow  into  or  out  of  the 
plate  as  functions  of  tlme^  the  system  Biot  numbers,  the  environment 
temperature  at  the  faces  of  the  plate,  and  the  initial  temperature  of 
the  plate . 

SECTION  III 
RESULTS 


3.1  GENERAL 

The  eigenvalues  were  calculated  from  Eq.  (12)  by  computer  for 
a  wide  range  of  characteristic  Biot  numbers  I'L  and  .  The  first  ten 

bj.  x52 

positive  roots  were  calculated  for  each  combination  of  s-nd  N_  , 

Bi  B2 

each  root  being  accurate  to  four  places.  Table  I  lists  values  of  e 
for  all  possible  combinations  of  the  following;  ^ 

W  =  0,  0.1,  0.2,  0.3,  0.4,  0.5,  1.0,  2.0,  4.0,  6.0,  8.0, 

^  10.0  ,  20.0,  100.0,  CO. 

W  =  0,  0.1,  0.2,  0.3,  0.4,  0.5,  1.0,  2.0,  4.0,  6.0,  8.0, 

^  10.0,  20.0,  100.0,  CO. 

The  series  in  Eq.  (ll)  and  (16)  converge  rapidly  with  ten  or 
fewer  roots  of  Eq.  (12)  for  all  values  of  Kt/i^  >  0.01. 

Unfortunately,  in  Eq.  (ll)  it  is  not  possible  to  obtain  only  the 
system  temperatiire  terms  on  the  left  side  of  the  equation  and  only  the 
Biot  and  Fourier  numbers  on  the  right  side  of  the  equation.  This 
eliminates  the  possibility  of  a  general  dimensionless  plot  of  Eq.  (ll). 
However,  if  it  is  assumed  that  the  initial  temperature  of  the  plate  V 
equals  the  environment  temperature  v^,  then  the  dimensionless  temperatture 
parameters  6  and  \[f  become  0  and  1,  respectively.  The  ri^t  side  of 
Eq.  (11)  then  becomes  a  function  of  N  ,  N  ,  and  Kt/f2  only.  This  is 

a  reasonable  assumption  for  some  of  the  typical  applications  that  were 
mentioned.  With  this  assumption,  the  quantitative  effect  of  varying 
conditions  of  surface  heat  transfer,  plate  thickness,  run  time,  and 
material  properties  on  the  temperatiire  response  of  the  plate  can  be 
determined.  It  should  be  emphasized  that  the  application  of  these 
equations  is  not  restricted  to  any  one  material,  environmental  tempera¬ 
ture  difference,  heat-transfer  coefficient,  or  run  time  because  of  the 
dimensionless  character  of  the  solution. 

3.2  ILLUSTRATIVE  EXAMPLE 

The  solution  can  best  be  appreciated  by  considering  a  typical 
problem.  Consider  the  case  of  an  exhaust  gas  diffuser  for  simulating 
the  hi^  altitudes  necessary  to  evaluate  the  perfcrmance  of, rocket  engines. 
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It  is  desired  to  know  the  tempera tiire  response  of  the  diffuser  wall 
hecause  the  response  must  he  limited  to  maintain  the  structural 
integrity  of  the  diffuser.  The  wall  is  initially  at  temperatin’e  V, 
cooled  hy  a  constant  temperature  =  V  water  reservoir  at  x  =  and 

at  t  >  0  is  subjected  to  a  high  temperature  gas  flow  at  temperature  v 

1 

at  X  =  0.  For  this  example  let  N_  be  0.3  and  0.6,  the  wall  be 

hi  hs  o 

of  3/S-in.  thick,  mild  steel,  the  cooling  water  temperature  be  70  r,  and 
the  hot  gas  temperature  be  4,000°F.  Determine  the  length  of  time  for 
the  wall  at  face  x  =  0  to  reach  800°F,  which  will  be  assumed  to  be  the 
limiting  structural  temperature. 

0  =  V  -  v^/v^-  =  1260-lih6o/530-Mi6o  =-32OO/-3930  =  0.8l4 


Referring  to  Fig.  2,  which  has  example  temperature  response  curves 
plotted  for  the  face  at  x  =  0,  gives  a  Fourier  number  Kt/X^  =0.4  for 

N_  =0.3  antd.  K,,  =0.6.  The  thermal  diffusivity  K  for  mild  steel  is 

0.49.  Substituting  gives 


0.49  t 


0.4 


0.4  (9.766  X  10“*)  =  7.972  X  lO"'*  =  0.000797  Hr 

0.49 


or  t  =  2.87  sec 

This  is  one  example  of  the  use  of  the  curves.  For  a  given  run 
time,  the  plate  temperature  could  have  been  determined  just  as  readily. 

For  other  given  conditions  the  plate  thickness,  coolant  water  temperature, 
or  coolant  side  Biot  number  can  be  found.  Figure  3  gives  example  surface 
temperatiu’e  response  curves  for  the  face  at  x  =  il.  From  Figs.  2  and  3, 
it  can  be  seen  that  the  magnitude  of  may  be  critical  in  reducing  the 

surface  temperature  response  of  the  wall.  For  example,  consider  the 
curve  =  0.2  in  Fig.  2.  For  a  given  material  wall  thickness 

and  run  time,  this  curve  represents  the  temperature  response  of  a  wall 
with  equal  heat-transfer  coefficients  B.t  tl/ Z  =  0  and  Vi/Z  =  1.  If  all 
conditions  remain  the  same  except  that  the  heat-transfer  coefficient  at 
x/i  =  1.0  is  quadrupled,  the  new  response  curve  has  the  Biot  numbers 
II_  =0.2  and  =0.8.  Depending  on  the  magnitude  of  the  Fourier 

number;  the  reduction  in  the  temperature  response  of  the  wall  may  or  may 

not  he  significant.  For  instance,  if  Kt/^^  =  1.0  the  difference  in 

the  response  of  the  curves  with  Biot  numbers  =  0.2  and 

hi  nz 

=  0.2,  W  =  0.8  is  about  four  percent;  however,  at  steady  state 

hi  132 

the  difference  is  about  35  percent.  Figure  3  shows  that  the  wall 
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temperature  at  x  =  responds  at  a  slower  rate  than  the  wall  temperature 
at  X  =  0^  as  one  would  expect.  Also,  the  temperature  response  decreases 
with  increasing  Ifg  for  a  fixed 

Figure  1+  gives  example  heat  storage  ctu:ves  for  various  combinations 

of  and  as  a  function  of  the  Fourier  number  Kt/i^. 

Hi  Hs 

SECTION  IV 
CONCLUSIONS 


Through  the  use  of  Eq.s.  (ll),  (l2),  and  (l6),  general  temperature¬ 
time  plots  and  energy-stored  plots  can  be  developed  to  cover  all  cases 
of  interest  for  a  given  situation.  For  design  purposes  the  temperature 
distribution  in  a  wall  is  of  importance  in  determining  thermal  stresses, 
structural  integrity,  and  peak  surface  temperatures .  For  the  example 

considered,  it  is  shown  that  the  ratio  of  K,,  to  can  be  of  significant 

bs  Hi 

importance  in  reducing  the  temperature  response  of  a  wall.  For  a 
specific  problem,  the  equations  may  be  used  to  determine  the  most  economical 
combination  of  wall  material,  wall  thickness,  and  coolant  flow  rate  and 
temperature j  or  even  if  it  is  feasible  to  limit  a  given  wall  to  an  acceptable 
temperature  response . 
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TABLE  I 

FIRST  TEN  POSITIVE  ROOTS  OF  THE  TRANSCENDENTAL  EQUATION  (fS  -  Nbi  Nb2)  TAN  («„)  =  (Nai  i  Nbj) 

FOR  VARIOUS  BIOT  NUMBER  COMBINATIONS 
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^3 

€4 

^5 

^6 

^7 
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.  C 

.1 
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3.1731 

6.2991 

9.4554 
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18.8549 
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.0 
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.0 

, 
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22.0093 
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•  5 
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.c 
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.  c 
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22.1697 
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.  G 

6 .0 
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4.1116 

6.9924 

9.9667 
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19.1531 

22.2545 
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28.4820 

.C 
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.0 
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,0 
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.  0 
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.  0 
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10.9956 
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26.7055 

29.6451 

.  1 

.  1 
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12.5823 
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.  1 

O 
•  C. 
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6. 3306 
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25. 
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.  1 

.3 

.6150 
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6.3462 
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22.0093 

25. 1486 
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.  1 

.4 
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3.2928 

6.3617 

9.4775 
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25.1526 
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.  1' 

.5 
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6. 3  77'1 
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2.0 
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.1 

4.0 
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.1 
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.2 
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.2 
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6.8551 
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,3 
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.3 
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.3 

10.0 

1.5979 

4.3690 
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16.2773 
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22.0546 
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.A 

2.0 
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18.9756 
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26.3588 
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12. 8981 
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22.3527 
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7.5461 
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22.0591 

25.1923 

28.3273 

.5 
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